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Abstract 

I discuss BBN with nonequilibrium v e <-> v s oscillations in the more general 
case of non-zero population of v s before oscillations 5N S ^ 0. I calculate 4 He 
primordial production Y P (5N S ) in the presence of v e <-> v s oscillations for different 
initial populations of the sterile neutrino state < 5N S < 1 and the full range of 
oscillation parameters. 

Non-zero 5N S has two-fold effect on 4 He: (i) it enhances the energy density and 
hence increases the cosmic expansion rate, leading to Y p overproduction and (ii) it 
suppresses the kinetic effects of oscillations on BBN, namely the effects on pre-BBN 
nucleon kinetics caused by the u e energy spectrum distortion and the v e — i7 e asym- 
metry generation by oscillations, leading to decreased Y p overproduction. Depending 
on oscillation parameters one or the other effect may dominate, causing correspond- 
ingly either a relaxation of the cosmological constraints or their strengthening with 
the increase of 5N S . 

I calculate more general BBN constraints on v e <-> u s oscillation parameters, 
corresponding to 3% Y p overproduction, for different initial populations of the sterile 
state. Previous BBN constraints were derived assuming empty sterile state before 
oscillations. The cosmological constraints for that case strengthen with the increase 
of 5N S value, the change being more considerable for nonresonant oscillations. 
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1 Introduction 



Big Bang Nucleosynthesis (BBN) provides one of the most sensitive probes of the physical 
conditions of the early Universe and is often used to constrain neutrino oscillations. 

Active-sterile neutrino oscillations may considerably effect the early Universe, and in 
particular BBN. They are capable of exciting additional light particles into equilibrium, 
thus, affecting the expansion rate H, and they may also distort neutrino energy spectrum 
and generate neutrino-antineutrino asymmetry, thus, in case of electron-sterile oscillations 
u e <-> u s , influencing the weak interaction rates of u e and, hence, the kinetics of nucleons 
at the pre-BBN epoch and the primordial nucleosynthesis. 

Among the light elements produced primordially, 4 He is the most abundantly pro- 
duced, most precisely measured and calculated [1, 2, 3, 4, 5], and usually it is the chosen 
element for putting BBN constraints on oscillations. Cosmological constraints on neutrino 
oscillations parameters, based on BBN and 4 He observational data, have been discussed 
in numerous publications, see for example refs. [6, 7, 8, 9, 10, 11, 12, 13, 14], the review 
papers [15, 16, 17] and the references there in. 

All available constraints were obtained assuming that the sterile state was empty 
before oscillations SN S = N u — 3 = 0, i.e. that it has been filled only due to active-sterile 
mixing during oscillations. N v is the number of neutrino species in equilibrium. 

In general, however, v s state may not be initially empty. There exist numerous pos- 
sibilities of producing v s before v e <-> v s oscillations: from GUT models, large extra 
dimensions models, Manyfold Universe models, mirror matter models, etc. They may 
be also produced in preceding T <-> v s oscillations in 4-neutrino and 5-neutrino mixing 
schemes. 2 The degree of initial population of u s , SN S , and its energy spectrum depends 
on the concrete model of u s production. 

Non-zero initially 5N S can influence the dynamical and the kinetic effects of oscillations 
on BBN [19] and change the cosmological limits on oscillations. The aim of this work 
is to define how and to what extend the available cosmological constraints on v e <-> v s 
oscillation parameters, obtained with the assumption 5N S = 0, are changed in the more 
general case of initially non-zero sterile population 5N S ^ 0. 

In the next section I discuss how 5N S ^ influences oscillations effects on BBN and 
provide qualitative considerations of the expected change of cosmological constraints. In 
section 3 I present the results of the numerical analysis of 4 He production in BBN with 

2 The analysis of experimental oscillations data provides some constraints on the sterile neutrino impact 
in the oscillations explaining atmospheric and solar neutrino anomalies [18], and hence, provides some 
indications about the possible values of SN S , which could have been eventually produced by oscillations. 
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nonequilibrium v e <-> v s oscillations and initially nonzero v s population < 5N S < 1 
and discuss the generalized cosmological constraints on oscillation parameters. In the 
last section the results and the assumptions of the model are resumed and possibilities of 
relaxing the cosmological constraints are discussed. 

2 Generalizing BBN Constraints on Oscillations 

I discuss the generalization of BBN constraints on v e <-> v s neutrino oscillations parameters 
for the case of oscillations effective after electron neutrino decoupling. It is known that at 
lower mass differences, when sterile neutrino production takes place after active neutrino 
decoupling, due to v e <-> v s oscillations between initially empty v s and electron neutrino, 
i.e. for (5m 2 /eV 2 ) sin 4 29 < 1CT 7 , the re-population of active neutrino becomes slow and 
kinetic equilibrium may be strongly broken by u e <-> v s oscillations. I will denote these 
oscillations with an initially non-equilibrium sterile neutrino further on 'nonequilibrium 
oscillations' for short. They may cause considerable deviations of the v e energy spectrum 
from the equilibrium Fermi-Dirac form [7, 10] (called further 'spectrum distortion') and 
generate neutrino- antineutrino asymmetry 3 , thus influencing the weak interaction rates 
of the processes governing pre-BBN nucleon kinetics u e +p <-> n + e + , u e + n^p + e,n^ 
p + e + u e the nucleons freezing and correspondingly the primordial elements production. 4 

For that oscillation case a numerical analysis of the evolution of the neutrino and 
nucleons, using the kinetic equations for the neutrino density matrix and neutrino number 
densities in momentum space was provided to make a proper account for the neutrino 
spectrum distortion, depletion and neutrino asymmetry growth at each momentum [10]. 
The analysis of such oscillations, accounting precisely for the kinetic effect, allowed to put 
stringent constraints on oscillations effective after active neutrino decoupling [10, 11, 16]. 

The analytical fits to the exact constraints, corresponding to 3% 4 He overproduction, 
are [10, 11, 16]: 

5m 2 (sin 2 2$) 4 < 1.5 x 10~ 9 eV 2 5m 2 > 
\5m 2 \ < 8.2 x 10~ 10 eV 2 5m 2 < 0, large 

3 The asymmetry effect at large mixings and small mass differences, discussed here, is subdominant [10] 
and leads to a slight suppression of the spectrum distortion effect. This is in contrast to the case of 
generation of lepton asymmetry by oscillations at relatively high dm 2 and small mixings [20], where it 
may play the dominant role. 

4 The effect of distortions on the energy distributions of neutrinos caused by residual interactions of 
neutrinos after 2 MeV during electron-positron annihilations was not considered, because it leads to a 
negligibly small change in Y p . [23] 
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These constraints, as well as all existing in literature BBN constraints, concern the 
case of empty sterile neutrino state at the start of oscillations. 

The presence of a nonzero SN S before oscillations exerts two types of effects on 
BBN [19]: 

(i) It increases the energy density by 5p = 7/8(T u /T^) 4 SN s p 1 , thus modifying the 
cosmic expansion rate H = ^j8irp/3Mp, which reflects into higher freezing temperature 
of the nucleons and overproduction of 4 He [21]. The dynamical effect of initially present 
SN S on primordial 4 He production will be denoted further on by 5Y d . Due to this effect 
strengthening of the cosmological bounds with respect to the ones calculated at 
SN S = should be expected. 

In case of u^ T <-> u s oscillations the dynamical effect is the only effect of non-zero 5N S 
present before oscillations. In the case of u e <-> v s oscillations with almost equilibrium 
neutrino energy distribution, i.e. oscillations taking place before neutrino decoupling, this 
is the leading effect as well. This effect can be accounted for simply by adding the initial 
5N S value to the one produced in oscillations. So, in both these cases the rescaling of the 
existing constraints is rather straightforward. 

(ii) In the nonequilibrium v e <-> v s oscillations case, the presence of partially populated 
v s suppresses the oscillations effects on pre-BBN nucleons kinetics [19]. Further on this 
kinetic effect is denoted by 5Y£, while the kinetic effects of oscillations are denoted by 
SYk (or in terms of the effective degrees of freedom 5N^, where 5Y^ ~ 0.013 x 5N^). 

Kinetic effects are a result of the generated energy spectrum distortion of v e in os- 
cillations between active and sterile neutrino, and to a smaller degree are due to the 
neutrino-antineutrino asymmetry, generated by oscillations. 5 Both are very sensitive to 
the degree of population of v s because the rate of oscillations is energy dependent. The 
spectrum distortion effect plays however the dominant role, so it will be discussed in more 
detail further on, although we have accounted precisely for both. 

The distortion of v e spectrum leads both to a depletion of the active neutrino number 
densities N u : N u ~ J dEE 2 n u (E) and a decrease of the r TO , causing an earlier n/p-freezing 
and an overproduction of ^e yield. The spectrum distortion is the greatest in the case 
the sterile state is empty at the start of oscillations, SN S = N Ue /N Us = 0. It decreases 
with the increase of the degree of population of the sterile state at the onset of oscillations 

5 Neutrino-antineutrino asymmetry generated during the resonant transfer of neutrinos exerts back 
effect on oscillating neutrino. Although its value is too small to have a direct kinetic effect on the 
synthesis of light elements, i.e. L « 0.01, it effects indirectly BBN by suppressing oscillations at 
small mixing angles, leading to less overproduction of He-4 compared to the case without the account of 
asymmetry growth [10]. 
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as illustrated in the following figure from ref. [19]. 



0.5 




Fig.l. The figure gives a snapshot of the spectrum distortion of the electron neutrino energy spectrum 
x 2 Pll(x), where x — E/T at a characteristic temperature 0.7 MeV, caused by resonant oscillations with 
a mass difference 5m 2 = 10~ 7 cV 2 and mixing sin 2 2d = O.f for different degrees of initial sterile neutrino 
population, namely SN S = (lower curve), SN S — 0.5 and SN S — 0.8 (upper curve). The dashed curve 
gives the equilibrium spectrum for comparison. 

So, larger SN S leads to a decrease of the kinetic effects of oscillations and to a decrease 
of the overproduction of 4 He by oscillations with respect to the case of initially zero 
SN S , i.e. <5iVfc < 5N® in , where SN® in is the kinetic effect corresponding to zero initial 
population of the sterile state, which presents in fact the maximal kinetic effect at a given 
set of oscillation parameters. 6 Due to the decrease 5Y£ in the overproduction of Y p , 
caused by the suppression of 5N k by the initially present SN S ^ 0, relaxation of the 
cosmological bounds with respect to the ones corresponding to an initially empty v s 
state should be expected. 

Obviously, there is an interplay between the two types of effects (i) and (ii) induced by 
non-zero 8N S . The total effect depends on the concrete values of the oscillation parameters 
and 5N S . The shift of the constraints may be expected in either direction (relaxing or 
constraining the existing constraints) and its analysis requires numerical study, accounting 
precisely for v e energy spectrum distortion and the suppression of the kinetic effects, which 
is essential. 

Using the approximate empirical formula giving the dependence of Y p production on 

6 For v e <-» v s oscillations effective after v e decoupling, these kinetic effects 5N% can be considerable, 
as large as 6N° = 6 [22]. 
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the two effects, obtained in ref. [19]: 



5Y P = 5Y d + 5Y k ° + 5Y k s ~ 0.013 x (5N S + 5N° km - 5N S x 5N° km ) 

it is possible to make some predictions concerning the value and the sign of the 4 He 
overproduction and the direction of the shift of the BBN. Here 5Y p = Y p — Y p stand , 5Yd = 
0.013 x 5N S , 5Y k ° = 0.013 x 5N% in , while 5Y k s = -0.013 x 5N S x 8N& n . 

First: The total effect of oscillations for initially non-zero SN S on 4 He pro- 
duction is smaller than the sum of the energy density increase effect (i) 5Yd 
and the maximum kinetic effect of oscillations (ii) 5Y k corresponding to zero SN S , 
due to the term SN S x SN kin expressing the decrease of oscillations kinetic effect with SN S . 

The results of the exact numerical study of the effects (i) and (ii) on 4 He overproduc- 
tion, illustrated in the following figures, confirm this estimation. In Fig. 2a and Fig. 2b the 
contribution of the different effects on neutron-to-nucleon freezing ratio X n = n^ n jn nuc is 
presented. (The Y p production mainly depends on X n , Y p ~ X n x exp(—t/r n ), where r n is 
the neutron lifetime.) The dotted curve shows the kinetic effect dependence on 5N S , the 
lower dashed curve gives the energy increase effect. The total effect is presented by the 
solid curve, which although has different behavior in the two cases (decreasing or increas- 
ing), is situated considerably lower than the uppermost long-dashed curve presenting the 
sum of the effects (i) and (ii). 
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Fig. 2a. The solid curve presents the frozen neutron number density relative to nuclcons X n — n£/n nttc 
as a function of the sterile neutrino initial population, at 5m = ±10~ 8 eV 2 , sin 2 29 = 1. The dotted 
curve presents the kinetic effect, while the lower dashed curve presents energy density increase effect. The 
uppermost long dashed curve corresponds to the total effect when the decrease of the kinetic effect is not 
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accounted for, as if the initial 5N S were stored in a state not participating in oscillations with electron 
neutrino, and the effects were simply additive. 

Hence, the cosmological constraints will be less stringent than in the case of simply 
additive effects. The latter case will have place when the enhancement of the energy 
density 5N S is due to other additional particles brought partially into equilibrium (like 
sterile neutrinos in the muon or tau-sectors, or other relic relativistic particles) while 
the sterile state, which participates further in oscillations with the electron neutrino is 
initially empty. In that case the two effects will be simply additive, and the uppermost 
long-dashed curve in Figs. 2 presents the overproduction of 4 He for such specific situation. 
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Fig. 2b. The solid curve presents the frozen neutron number density relative to nucleons — ii{/u„„ c 
as a function of the sterile neutrino initial population, at 5m = ±10~ 9 eV 2 , sin 2 29 = 1. The dotted 
curve presents the kinetic effect, while the lower dashed curve presents energy density increase effect. 
The uppermost long dashed curve corresponds to the total effect if the effects were simply additive. 

However, due to the fact that SN kin is a decreasing function of 5N S , naively adding 
the two effects exaggerates Y p overproduction, and hence would define stronger bounds 
than the real ones. 

Second: The direction of the shift of the constraints should be as follows: 
in case 5N® in 5N s < 5N S the constraints will be strengthened in comparison with 
the 5N S = constraints, while in the opposite case they will be relaxed. 

In fig. 2a the solid curve, presenting the total effect, is a decreasing function of 8N S 
because 5Njt in > 1, i.e. for that set of oscillation parameters the overproduction of 4 He due 
to oscillations decreases with the increase of the initial population of the sterile neutrino 
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state. Obviously the suppression effect (ii) dominates. In that case we expect relaxation 
of the cosmological constraints compared to the case of initially zero population of v s . 

Fig. 2b presents the results for a set of oscillation parameters for which 5N® in < 1, then 
as it is illustrated the total effect is an increasing function of 5N S , the dynamical effect 
(i) dominates over (ii), thus the overproduction of 4 He increases in respect to the case of 
5N S = and hence, the cosmological constraints must become more stringent with the 
increase of 5N S . And as far as 5N® in = SN tot at SN S = 0, the BBN constraints, corre- 
sponding to 4 He observational uncertainty expressed as SN tot < 1, will be strengthened. 
For example, 4 He uncertainty 5Y P ~ 0.007 corresponds to SN tot ~ 0.54. Thus at 5N S = 
5Nj* in ~ 0.54 < 1, and hence strengthening of the cosmological constraints is expected. 

The exact form of the cosmological constraints obtained by a detail numerical study 
of the effects (i) and (ii) and corresponding to such 4 He overproduction and different 5N S 
values is presented in the next section. 

In case of bigger 4 He uncertainty, parametrized by 5N tot > 1, the term (ii) dominates 
(as illustrated in Fig. 2a), leading to a decrease of the 4 He overproduction due to oscil- 
lations and hence, relaxation of the cosmological constraints should be expected. 7 It is 
interesting to note that contrary to some prejudice, even for 4 He uncertainty equivalent 
to SN tot > 1 cosmological constraints on oscillation parameters still persist, provided that 
5 N s < 1 and that a proper description of the neutrino energy spectrum distortion is made. 

3 BBN constraints for partially filled v s 

For calculating cosmological constraints on oscillation parameters, I have followed the 
kinetic approach described in detail in ref. [19]. The case of v e <-> v s oscillations effective 
after electron neutrino decoupling is considered. For simplicity mixing just in the electron 
sector is assumed. The more general case of a sterile neutrino state being partially filled 
before the oscillations become effective is studied. In a previous work [19] it was shown 
that for a wide range of 5N S values the kinetic effects of oscillations, namely spectrum 
distortion of the electron neutrino and lepton asymmetry growth, play a considerable 
role, hence rough analytical estimations of the kinetic effects are not applicable [25, 10]. 
Therefore, here an exact numerical analysis of 5N S effects (i) and (ii) on pre-BBN nucleons 
freezing and the production of 4 He in BBN with oscillations is provided. 

The kinetic equations for the density matrix of neutrino and antineutrino ensembles 
are simultaneously and selfconsistently solved with the kinetic equations governing the 

7 In more detail this case will be discussed in a following paper [24] . 



7 



evolutions of nucleons from the period of neutrino decoupling T ~ 2 MeV till the nucleons 
freezing epoch, down to T ~ 0.3 MeV. The neutrino spectrum is described by 1000 bins 
for the nonresonant case, and by more than 5000 bins (depending on the asymmetry 
growth region of oscillation parameters). This detail study, performed for a wider range 
of parameters, confirms the importance of the kinetic effect at non-zero 5N S , found in 
ref. [19]. I calculate primordial 4 He yield Y p (5N s , 5m 2 , sin 2 2d) at different 5N S values for 
the full set of oscillations parameters of the model: for all mixing angles $ and mass 
differences 5m 2 < 10~ 7 eV 2 . 

For the analysis of the constraints the observational uncertainty of the primordially 
produced 4 He is assumed 5Y P < 0.007 in correspondence with the accepted systematic 
error in the Tie measurements [1]. Then the maximum possible value of SN S at BBN 
epoch is constrained on the basis of BBN considerations: Using the approximate empirical 
formula SY p ~ 0.0138N tot: 5Y P < 0.007 corresponds to 5N S < 0.54. So, in our analysis 
we have varied SN S in the range 0.0 < SN S < 0.5 with a step 0.1. The case SN S > 0.54 
corresponds to higher "tie uncertainty and will not be studied here. 

This choice of maximum 8N S is also supported by the results of the recent standard 
BBN analysis + Tie observational data with an input the baryon density value from 
WMAP data, which provide bounds on the number of additional neutrino species at 
BBN in the range 5N S < 0.1 - 0.5 [2, 3, 26, 27, 28]. 

Uncertainty 8Y p = 0.007 corresponds to 5Y P /Y P ~ 3%, so in this analysis, I have 
assumed that primordial 4 He abundance is known with an accuracy better than 3%, and 
will call the isohelium contours corresponding to 5Y P = 0.007 '3% 4 He overproduction 
contours', and the obtained BBN constraints '3% 4 He constraints'. 

In Fig. 2 I present the calculated cosmological constraints on oscillation parameters 
for the 5N S = 0.1, and 5N S = 0.5 for the resonant (to the left) and the non-resonant 
(to the right) oscillation cases. The region upwards of the corresponding curves is the 
cosmologically excluded. The upper curve is for 5N S = 0.1 case, the lower one - for 
8N S = 0.5. 5N S = constraints are given for comparison by the dashed contours. 
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Figure 3: BBN constraints on oscillation parameters for the resonant (l.h.s.) and the non- 
resonant v e <-> u s oscillations and for initial degrees of population of the sterile neutrino state 
SN S = 0.1 and SN S = 0.5. The dashed contours present the constraints for 5N S = case for 
comparison. 

The analysis shows that the two effects (i) and (ii) of 5N S , nearly compensate each 
other for small 5N S values. Hence, the cosmological constraints for SN S = 0.1 slightly differ 
from the ones for 5N S = 0, as illustrated in Fig. 3. For SN S > 0.1 SYj dominates over the 
suppression term 5Y£, so the constraints are strengthened. As a whole the cosmological 
constraints in the resonant case are slightly changed compared to the case with initially 
zero sterile state, while in the non-resonant case the change is more noticeable, and the 
constraints are becoming more stringent with the increase of SN S . 

The precise account of the kinetic effects of oscillations on BBN, i.e. the study of 
the distorted spectrum distribution of neutrinos and its effect on nucleons kinetics, allows 
strengthening the cosmological constraints by an order of magnitude towards smaller mass 
differences in comparison with calculations considering just the integral effect like neutrino 
number density depletion. 
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4 Discussion 



The presence of a non-empty sterile state before v e <-> v s oscillations was not considered 
in previous cosmological constraints on oscillation parameters. This work is a step to- 
wards generalizing the cosmological constraints on oscillation parameters of neutrino. It 
discusses the effect of partially filled sterile state before oscillations on primordial produc- 
tion of 4 He and on the BBN oscillations constraints. 

It is shown that initially non-zero u s state has two fold effect on 4 He: on one hand 
it affects Universe dynamics leading to overproduction of Y p , while on the other hand its 
presence suppresses oscillations kinetic effects causing a decrease of Y p overproduction. 
Thus, depending which effect dominates, non-zero initially v s may cause strengthening or 
relaxation of the cosmological constraints on oscillation parameters. 

I have calculated the isohelium contours corresponding to 3% overproduction of "fte, 
in a model of BBN with electron-sterile oscillations effective after neutrino decoupling 
and for the general case of the sterile neutrino state being partially populated before 
oscillations, 0.0 < 5N S < 0.54. The constraints become more stringent with the increase 
of SN S value. So, the presence of non-zero initial population 5N S < 0.54 strengthens the 
3% ^e cosmological constraints. 

These cosmological constraints are obtained for the case of 2-neutrino mixing. Further 
generalization of the constraints includes the account of mixing between active neutrinos, 
which has been proved important in the resonant oscillation case for oscillations proceed- 
ing before neutrino decoupling [14]. 

The constraints are obtained for the natural assumption that at BBN epoch the initial 
lepton asymmetry is of the order of the baryon one. However, in fact the lepton asymmetry 
in the neutrino sector is not strongly constrained [29]. We expect that, provided a small 
lepton asymmetry L « 0.1 is present, it may be large enough to relax or alleviate 
the discussed BBN bounds on neutrino oscillations. For 5N S = case it was proven 
that larger than 10~ 7 lepton asymmetry may strongly suppress oscillations effective after 
neutrino decoupling and change the constraints, while initial lepton asymmetry larger 
than 10~ 5 10" 4 can alleviate them [30, 11]. 

Cosmological constraints can be relaxed also if the systematic error of Y p is higher 
than 0.007. However, it is interesting to note that even for SY p ~ 0.01, and a considerable 
initial population of u s 5N S < 1, the constraints may not be removed, but just relaxed, 
as the results of ref. [24] suggest. 
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